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ABSTRACT 

We present axisymmetric, orbit-based models to study the central black hole, stellar mass-to-light 
ratio, and dark matter halo of NGC 4594 (M104, the Sombrero Galaxy). For stellar kinematics, we 
use published high-resolution kinematics of the central region taken with the Hubble Space Telescope, 
newly obtained Gemini long-slit spectra of the major axis, and integral field kinematics from the 
SAURON instrument. At large radii, we use globular cluster kinematics to trace the mass profile and 
apply extra leverage to recovering the dark matter halo parameters. We find a black hole of mass 
M,= (6.6 ± 0.4) X 10^ Mq, and determine the stellar M/Li = 3.4 ± 0.05 (uncertainties are the 68% 
confidence band marginalized over the other parameters). Our best fit dark matter halo is a cored 
logarithmic model with asymptotic circular speed Vc — 376±12 km s ^^ and core radius r^ = 4.7±0.6 
kpc. The fraction of dark to total mass contained within the half-light radius is 0.52. Taking the 
bulge and disk components into account in our calculation of (Je puts NGC 4594 squarely on the M-a 
relation. We also determine that NGC 4594 lies directly on the M-L relation. 

Subject headings: galaxies: individual (M104, NGC 4594) — galaxies: kinematics and dynamics — 
galaxies: photometry 



1. INTRODUCTION 

Most galaxies are thought to host supermas- 
sive black holes (SMBHs) at their centers. The 
masses of these SMBHs have been observed to cor- 
relate with several properties of their host ellip- 
tical galaxies and of the classical bulge compo- 
nents of their host disk galaxies. For example, 
M, correlates wi th ga l axy/bulge mass (iDress lcr 1 98ft 
Magorrian et all 119981: iLaor 2001: iMcLure fc Dunlojol 
2002: Marconi fc HuntI 12003 : Haring fc Rixl 12004 , 
luminosity (the M-L relat i on) (Kormcndy 199^ 
Kormen dv fc Richstond 119951 : iKormendv fc Gebhardti 



200lH Giiltekin et al.ll2009bD. velocity di s persion (the M- 
a relation) (Ferrarcse & Mcrritt 2000: 'Gcbhardt et al.l 
[2000; Trcmainc et al. 2002; Giiltckin et al. 2009 ^, and 
globular cluster con tent ([Burkert fc Tremaind 120101 : 
[Harris fc Harrid 120111) . These and other, similar cor- 
relations suggest that galaxy formation and black hole 
growth are fundamentally linked. To better understand 
this interplay, accurate black hole masses are needed. 

One challenge that limits the accuracy is the deter- 
mination of the host galaxy's inclination. Projection ef- 
fects are difficult to model and cause loss of information, 
leading to systematic uncertainties. Therefore, SMBHs 
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in galaxies whose inclination is confidently known have 
the best chance of being accurately and robustly mea- 
sured. Another issue that limits the accuracy is the 
effect a dark matt er halo has on th e deter mination of 
SMBH mass. Ge bhardt fc Thomad (|2009f ) show that 
orbit-based models can underestimate black hole mass 
whe n dark matter is not consider ed in the modeling, how- 
ever iSchulze fc Gebhardtl (|2011| ) find the effect is small 
when the black hole's sphere of influence is well-resolved. 
NGC 4594 (M104, or the Sombrero Galaxy) is a 
nearly edge-on Sa type spiral with a prominent stellar 
disk and large, classical bulge (Kormcndy & Kcnnicutl 
120041) . The shape of this disk indicates that it (and 
thus the entire galaxy) is inclined at an angle very 
close to 90°. Throughout this paper we assume a dis- 
tance to NGC 4594 of 9.8 M pc, calculated f rom sur- 
face brightness fluctuations ( Tonrv et al.ll200"ll ). Unless 
otherwise stated, all distance-dependen t quantities are 
scaled to this value. iTonrv et al.l ()2001f ) use a value of 
Hq = 74 km s^ Mpc~^ in their distance determinations 



however we compare our M, and Ly to iGiiltekin et al.l 
(2009b) who adopt i?n = 70 in the ir work. We there- 
fore scale the IGiiltekin et all ()2009bD distances down by 
6%. Black hole mass scales as M.oc D and luminosity as 
Ly oc D~^; these quantities are adjusted accordingly. 

NGC 4594 was one of the first galaxies in which a black 
hole was d etected, and it has a long history of study. 
lKormendvl (1988. hereafter lK88i ) first found evidence for 
a massive black hole of M,= 5A^]^\^ x 10^ M© using 
only gr ound-based observation s. With isotropic Jeans 
models. [Emsellem et al.l ()1994b[ ) measured a black hole of 
mass M,= (5.4± 0.5) x 10* Mq. Later. IKormendv et all 
(1996, hereafter 'K96) used high-resolution kinematics 
from the Faint Object Spectrograph (EOS) on the Hubble 
Space Telescope (HST) — the same data set we include — 
to measure log M, = 8.8 ± 0.5 Mq. This corresponds 



to a mass of 5.8l4^Q x IO^A/q. With isotropic mod- 
els, iMagorrian et al.l ()1998l ) obtained a value of A/,= 
6.9^0.? X 10^ Mq. These values for M, aU lie towards the 
high mass end of the M-a and M-L relations. Massive 
SMBH measurements are frequently being revised, and 
we expect the confidently known inclination of NGC 4594 
to lead to one of the more secure measurements of a high 
mass SMBH. 

We present new Gemini spectroscopy of the major axis, 
as well as SAURON integral field kinematics covering the 
central region of the galaxy. We also use high-resolution 
HST /FOS kinematics of the nucleus and kinematics de- 
rived from globular clusters at large radii. We com- 
bine these kinematic datasets with HST and ground- 
based photometry to run axisymmetric orbit-based mod- 
els. These models allow us to measure the black hole 
mass, stellar mass-to-light ratio, and dark matter halo 
of NGC 4594. In addition, we also recover information 
about the internal orbit structure of the galaxy. 

2. DATA REDUCTION AND ANALYSIS 

Dynamical modeling requires as input the three- 
dimensional luminosity density distribution ^{r), as well 
as the line-of-sight velocity distribution (LOSVD) at 
many locations in the galaxy. We use HST and ground- 
based images for the photometry. Our kinematics include 
high-resolution HST /FOS spectra, long-slit spectra from 
GNIRS on Gemini, SAURON integral field kinematics, 
and individual velocities of globular clusters (GCs). We 
discuss each in turn. 

2.1. Photometry 

In order to cover a large enough dynamical range to 
have leverage on both the central SMBH and dark halo, 
we use surface brightness profiles from HST and ground- 
based images. The stellar disk of NGC 4594 domi- 
nates at intermediate radii on the major axis causing 
the isophotes in this region to be substantially fiattened. 
This abrupt change in ellipticity introduces an additional 
challenge to the deprojection. Our standard technique 
is to assume that the surfaces of cons tant luminosity 
density i/ are coaxial, similar spheroids (jGebhardt et al.l 
119961) . Clearly the presence of a disk invalidates this as- 
sumption, so we decompose the surface brightness into 
bulge and disk components, deprojecting each separately 
so that our assumption holds for each component. Af- 
terwards, we re-combine the deprojected profile of each 
component t^uigc + ^^disk and input the total z^(r) into our 
modeling program. 

The bulge-disk decomposition fits directly to a pro- 
j ected image. We construct a model disk by considering 
a lSersid(|1968[ ) profile: 



fi{R) = ^oexp[-(i?/i?o 



,1/r; 



(1) 



where /io is the surface brightness at R — Rq and n is 
the Sersic index. For n = 1, the profile is an exponential. 
For inclinations other than 90°, the projection of our 
disk model is an ellipse. By specifying the inclination of 
the disk i, the axial ratio b/a of the ellipse is given by 
b/a = cos i for a thin disk. 

We construct many disk models by varying /ig, ^O; 
i, and n (keeping n close to 1). Each model is then 
subtracted from the image until the residual brightness 



TABLE 1 
Summary of Disk Parameters Fit 






Disk fio (mag arcsec"^) Ro (arcsec) 


n 


i 


Outer 1 18.8 66.8 
Outer 2 16.7 40.1 
Nuclear 20.4 4.1 


1.0 
1.0 
1.1 


80 
80 
83 



distribution has elliptical isophotes. The remaining light 
is assigned to the bulge component. A ID major axis 
bulge profile is produced by averaging the bulge light in 
elliptical, annular isophotes. Hence, we are left with an 
analytic disk model and a non-parametric bulge model. 
We identify the best bulge and disk models as those that 
minimize the rms residuals of the model-subtracted im- 
age. 

In addition to the obvious main disk, NGC 4594 hosts 
a well-studied nuclear disk (Burkhcad 1986,K88,K96) at 
small radii. We fit the nuclear disk in the HST image 
and the main disk in the ground-based image. Because 
we fit directly to the images, dust lanes and object mask- 
ing become important. We keep a bad pixel list which 
instructs our code to ignore trouble spots. Dust lanes 
are s elected by eye, while SExtractor (jBertin fc Arnouti 
Il996f ) is used to identify foreground stars, background 
galaxies, and globular clusters. 

2.1.1. HST Image 

To probe the nuclear region, we use a PSF-deconvolved 
HST Wide Field Planetary Camera 2 (WFPC2) i mage 
(GO-5512; PI: Faber). This image is presented in IK96I 
and provides an excellent view of the central region of the 
galaxy. Centered on the PCI camera, the image is taken 
in the F547M fiher, and has a scale of 0'.'0455 pixel"^ of 
the central 34" x 34" of the galaxy. The P SF deconvo- 
lution uses the L ucy- Richardson algorithm (|Richardsonl 
Il972t ILucvlll97^ for 40 iterations and is well-tested on 
WFPC2 images (jLauer et al.lll998D . The best fit parame- 
ters from the bulge-disk decomposition are listed in Table 

m 

NGC 45 94 is als o thought to have weak LINER emis- 
sion (Bendo e t al.l [2006) and there is a point source in 
the HST image. Furthermore, heavy dust absorption 
also makes the determination of the central bulge surface 
brightness profile difficult for R < C'l?. To deal with 
these issues, we extrapolate the bulge surface brightness 
/^buigc(-R) inwards to R — G"02 with a constant slope fit 
to the region near R — 0"17. Figure [1] shows the result 
of this extrapolation, as well as the other components fit 
in the ground-based image. 

2.1.2. Ground-Based Image 

We obtained a wide-field, /-band image from the Prime 
Focus Camera on the McDonald 0.8 m telescope. This 
instrument provides a large unvignetted field of view 
(45x45 arcmin^ ) and a single CCD detector. Therefore, 
we can more robustly carry out sky-subtraction and ac- 
curately constrain the faint isophotes. The image is cor- 
rected for bias, fiat field, and illumination using standard 
routines in IRAF. 

In our fit to the ground-based image, we ignore the 
central 20" due to over-exposure and contamination from 
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Fig. 1. — Radial plot of all resulting components from our bulge- 
disk decomposition along the major axis. Dotted lines are the disk 
profiles with parameters from Table [T] The solid black line is the 
total surface brightness (bulge + disks). Solid colored lines are the 
bulge profiles, the red line is the result from fits to the HST im- 
age, and the blue line is from the ground-based image. Diamonds 
indicate the raw HST bulge profile before we apply our dust cor- 
rection and point-source removal. The gap between the HST and 
ground-based bulge profiles is interpolated over before deprojec- 
tion. Plotted in green is the globular cluster surface brightness 
profile, arbitrarily scaled to match the stellar surface brightness at 
its innermost point. 

the nuclear disk. We attempt fits in the region 20" — 900" 
with only one stellar disk given by Equation ^ , however 
these produce unacceptable residuals. Instead of modify- 
ing Equation ([T]) , we add a second disk (in addition to the 
nuclear disk fit only in the HST image). This approach is 
similar to the Multi- Gaussian Expansion technique used 
to model the light dist ribution of bulges and ellipticals 
(jEmsellem et al.lll994af) . A summary of all components 
fit is given in Table [T] and plotted in Figure [1] 

2.1.3. Globular Cluster Profile 

Globular clusters (GCs) are essentially bright test par- 
ticles that allow us to probe the potential at radii where 
the stellar light is faint. They ha ve been used in orbit- 
based models of other galaxies (iGebhardt fc Thomas! 
[20091: IShaTfe Gebhardtl[20T0l: IMurphv et al.ll201lD . To 
include the m in our models, we use the GC number den- 
sity profile (jRhode fc Zepj[20M ) as an analog to the stel- 
lar density. The number density profile is converted to 
a surface brightness profile by arbitrarily adjusting the 
zero point to match the stellar profile in log space. 

The green line in Figure [T] shows that the slope of the 
GC surface brightness profile is different from that of the 
stars. We run models using both the measured luminos- 
ity density distribution of the GCs and assuming that of 
the stars. We find significant preference for the measured 
GC profile. 

2.1.4. Bulge Profile 

Our bulge-disk decomposition returns a non- 
parametric form of the bulge profile. It is not 
necessary to have a parameterized bulge profile for our 
dynamical models, however we fit a Sersic profile to 
ground-based bulge model using Equation ([1]). The 
bulge is well-fit by a Sersic function, with the rms 
residuals equal to 0.08 mag arcsec^^. We measure 



^0 = 13.5 mag arcsec~^,i?o — O'.'l, and n = 3.7. 

We can convert the central surface brightness /io and 
radius Rq parameters to the more familiar "effective" 
parameters /ie and Re. The effective radius Re is given 
by i?e = {bn)^Ra and the eff ective surfac e brightness 
fie ^ Ho + 2.51og(e)&„ (jMac Arthur et all [2003) . The 
fa ctor bn depends on n; an expansion for bn can be found 
in lMacArthur et"aLl ()2003[ ). Applying these conversions, 
we obtain fXe = 21.3 mag arcsec"^ and Re — 156.2". 

We obtain a simpler estimate for the half-light radius 
of the bulge from integration of the surface brightness 
profile; no fitting functions are required. We estimate 
Re = 117" ± 12". The integrated magnitudes are cal- 
culated for component x by L^ = 2Tr J Ix{r)r dr. This 
does not take into account the ellipticity of each compo- 
nent, so we scale the luminosity by L^J"^^ « (1 — ex)Lx 
where tx is the ellipticity of each component, assumed 
to be constant with radius. The bulge profile is known 
to become rapidly circular for r > 100" (jBurkheadlll986[) 
so our procedure almost certainly underestimates Mbuigo 
and B /T . These numbers are computed as a sanity check 
only, and do not affect the dynamical models. 

The absolute integrated magnitudes are Mdisk = —21.4 
and Mbuige = —22.5 in F8 1 4 W, corrected for Galact ic ex- 
tinction along the line of sight (jSchlegel et al.llT998f) . Us- 
ing the HST cahbration package SYNPHOT (described 
in detail below), we convert these F814W magnitudes to 
y-band Vega magnitudes. We obtain My^buigc = —22.1 
and My^disk — —21.0. These structural parameters lie 
exactly on the fundame ntal plane for b u lges and el- 
lipticals as presented in IKormendv et al.l (j2009() . Our 
integrated magnitudes translate to a bulge-to-total ra- 
tio B/T = 0.73 with the nuclear disk contributing 1% 
of the total light. Thi s value of B/T i s low er than 
previous measurements — IKormendv et aLl ()2011[ ) report 
B/T = 0.925 ± 0.013. Our B/T, howe ver, is in good 
agreement with a recent measurement by IGadottil (J20111 
model ED). Regardless of the value of B/T, our dynam- 
ical models are unaffected, because we add all the bulge 
and disk light together again after the deprojection. 

We do not explore the possibility of fitting an expo- 
nential stellar halo in addition to a bulge and disk as 
IGadottil ()2011D do. Our bulge-disk decomposition pro- 
duces a non-parametric bulge profile which could in prin- 
ciple be a combination of a Sersic bulge plus exponential 
halo. However, this resulting profile is well-fit by a Sersic 
function with n significantly larger t han 1. We ther efore 
do not agree with the claim made bv lGadottl ()2011[ ) that 
the bulge of NGC 4594 is actually an exponential stellar 
halo. 

2.1.5. Deprojection 

We combine the HST and ground-based bulge pro- 
files by zero-pointing both to F814W. We calculate the 
F547M photometric zero point for the HST image from 
the SYNPHO T pac kage in IRAF. Spectral template fit- 
ting (Section 12.2.11) shows that in the central region of 
the galaxy > 85% of the light comes from K6III stars. 
We therefore convert the F547M zero point to F814W 
with SYNPHOT using the Bruzual AtlaE0 template for 
a K6III star. 

^ http://www.stsci.edu/hst/observatory/cdbs/bz77.html 
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Fig. 2. — Example Gemini spectrum. Top: observed galaxy spec- 
trum (black) and best-fit LOSVD-convolved template star (red). 
Dotted lines indicate regions of the spectrum ignored in the fit. 
Bottom: spectrum of the template star. The velocity dispersion of 
the LOSVD in this fit is cr = 190 ± 12km s'^ 

Before deprojection, we extrapolate the one- 
dimensional profiles /i(i?) with a constant slope to 
R = 1800". The three disk profiles are then combined 
and deproje cted via Abel i nversio n in the manner 
described in iGebhardt et alJ ()1996[ ). We assume an 
inclination of i = 90°. The inclinations of the combined 
disk components imply an ellipticity of e = 0.83. Our 
composite bulge profile is deprojected in a si milar fash- 
ion, assuming a constant ellipticity of 0.25 (jBurkheadI 
Il986( ). We then add fdisk('', ^) + t'buigG(?', ^) to obtain 
the total luminosity density distribution iy(r, 0) input to 
our models. 

The globular cluster luminosity density profile is ob- 
tained via a similar deprojection, but with the additional 
assumption of spherical symmetry. The normalization of 
the GC light profile is irrelevant, as our models fit only 
to the slope of the profile. 

2.2. Kinematics 

Kinematics for NGC 4594 come from four sources. The 
first uses near-IR data from Gemini/ GNIRS long-slit ob- 
servations along the major axis. These data were taken 
under good seeing conditions (around 0'.'5) and has high 
S/N. The second set comes from the Faint Object Spec- 
trograph (FOS) on HST using the s quare aperture of 
0'.'21 X 0'.'21 and is pubhshed in lK96i The third set of 
data is from the SAURON instrument (jEmsellem et al.l 
[2003). The SAURON data for NGC 4594 have not been 
published previously. Individual velocities from globular 



clusters are our fourth source of kine matics. These data 
are published in [Bridges et al.l (|2007[ ). We describe each 
dataset in detail. 

2.2.1. Gemini Kinematics 



We use GNIRS (jElias et al.ll2T)0fil ) on the Gemini South 
Telescope to measure near-IR spectra of NGC 4594. The 
data were taken on 17 January, 2005. We placed the 
150"x0'.'30 slit along the major axis with the galaxy nu- 
cleus centered within the slit. We use a spatial pixel size 
of 0'.'15. With the 321 1/mm grating in 3rd order, we ob- 
tain a wavelength coverage of 19800-26200 A at 6.4 A 
per pixel. Using sky lines, we measure a resolving power 
around 1700 or an instrumental dispersion of 75 km s^^. 
The total on-target exposure is 24 minutes, taken in 12, 
2-minutc individual exposures. Sky frames of equal ex- 
posure are taken throughout. 

From both setup images and images of telluric stan- 
dards, we measure a FWHM in the spatial direction of 
0'.'5, assuming a Gaussian distribution. We use this PSF 
for the dynamical models. 

We use a custom pipeline to reduce the GNIRS data; 
however, it produces very similar results to the Gemini 
GNIRS reduction package. The pipeline includes dark 
subtraction, wavelength calibration for the individual ex- 
posures, sky subtraction, registration and summing. 

There is adequate signal to extract kinematics out to 
a radius of 45". Figure [2] shows an example spectrum, 
where we plot the data in black and the template con- 
volved with the best-fit line-of-sight velocity distribution 
(LOSVD). The velocity templates come from the GNIRS 
spectral library (jWinge et al.f 2009') , where we select stars 
with a range of types from G dwarf to late giant. The 
kinematic extraction program performs a simultaneous 
fit to the LOSVD and relative w eights of the templates . 
Thi s procedure is descri bed in IGebhardt et al.l ()2000l ) 
and iPinknev et al.l (120031 ) . We present these data in the 
form of Gauss- Hermite moments in Table [21 

Figure [3] shows the kinematics derived from our anal- 
ysis of the Gemini spectra. Between 1" and 5", V rises 
and a drops. This is the result of the nuc lear d isk which 
becomes important at this radial range (jK88l ). Beyond 
10", we see similar behavior in V and a, it is caused by 
the main stellar disk. 

2.2.2. HST/FOS Kinematics 

IK96I present HST/FOS kinematics of the nuclear re- 
gion of NGC 4594. The FOS has a 0'.'21 x 0'.'21 aperture. 
There are three pointings with accurately known posi- 
tions for NGC 4594 (GO-5512; PI: Faber). The dynam- 
ical models include the e xact placement and aperture 
size for the FOS point ing (|K9l ). and use the HST PSF 
(IGebhardt et al.l[2000l) . 

2.2.3. SAURON Kinematics 

We also include SAURON integral field kinematics The 
SAURON data are from a single pointing exposing on the 
central regi on, taken in the low resolution setting of the 
instrument (jBacon et al.ll2001l ). In addition to V and a, 
the SAURON data also include the higher order Gauss- 
Hermite moments hs and /i4. Detai ls of the data reduc - 
tion and a nalysis can be found in I Bacon et al.l ()2001l) ; 
lEmsellem et al- (,200l ). 
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TABLE 2 
Gemini Kinematics 



R arcsec 


V kms-i 


AV km s-i 


(T km s ^ 


Act km s~^ 


^3 


A/l3 


/l4 


Ahi 


0.00 


19 


14 


253 


16 


-0.087 


0.033 


0.023 


0.046 


0.15 


-36 


11 


257 


12 


-0.008 


0.042 


-0.016 


0.038 


0.30 


-75 


12 


249 


8 


-0.053 


0.048 


-0.017 


0.037 


0.52 


-112 


11 


234 


8 


0.020 


0.037 


-0.047 


0.035 


0.82 


-144 


12 


221 


9 


0.060 


0.038 


-0.051 


0.032 


1.20 


-172 


9 


202 


9 


0.065 


0.039 


0.003 


0.031 


1.73 


-190 


8 


185 


9 


0.107 


0.031 


0.018 


0.031 


2.40 


-208 


7 


184 


9 


0.089 


0.032 


0.024 


0.031 


3.30 


-232 


8 


175 


9 


0.140 


0.029 


0.044 


0.027 


4.57 


-236 


9 


171 


7 


0.173 


0.032 


0.023 


0.027 


6.45 


-235 


7 


178 


9 


0.165 


0.033 


0.058 


0.025 


8.77 


-189 


11 


203 


11 


0.028 


0.039 


-0.007 


0.036 


11.40 


-171 


10 


192 


10 


0.071 


0.033 


-0.009 


0.033 


14.32 


-187 


16 


185 


15 


0.227 


0.042 


0.102 


0.060 


17.70 


-201 


13 


229 


16 


-0.023 


0.042 


0.020 


0.050 


22.20 


-228 


12 


198 


14 


0.125 


0.040 


0.052 


0.047 


28.58 


-235 


11 


194 


14 


0.122 


0.042 


0.036 


0.051 


36.08 


-285 


6 


141 


10 


0.046 


0.041 


0.032 


0.034 


44.40 


-277 


7 


149 


9 


0.094 


0.042 


0.101 


0.039 


-0.15 


45 


12 


240 


15 


-0.026 


0.034 


0.002 


0.041 


-0.30 


112 


11 


243 


15 


-0.052 


0.036 


-0.016 


0.038 


-0.52 


130 


10 


224 


15 


-0.080 


0.043 


0.014 


0.037 


-0.82 


162 


8 


212 


10 


-0.068 


0.049 


0.019 


0.030 


-1.20 


176 


7 


206 


11 


-0.066 


0.047 


0.065 


0.031 


-1.73 


205 


7 


190 


12 


-0.087 


0.043 


0.022 


0.037 


-2.40 


226 


7 


173 


12 


-0.101 


0.037 


0.053 


0.039 


-3.30 


244 


9 


187 


11 


-0.113 


0.039 


0.037 


0.037 


-4.65 


246 


8 


184 


12 


-0.130 
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Note. — Kinematics along the major axis of NGC 4594. Gauss-Hermite moments were derived 
from the LOSVDs that are the input to the dynamical models. 



Our modeling code fits to the entire LOSVD rather 
than its moments, so we reconstruct LOSVDs from the 
Gauss-Hermite moments. We create 100 Monte Carlo 
realizations of a non-parametric LOSVD from the un- 
certainties in t he Gauss-Hermite parame ters of each 
SAURON bin (|Gebhardt fc ThomasI 120091) . The 1433 
reconstructed SAURON LOSVDs are spatially sampled 
more finely than our modeling bins. We therefore average 
the SAURON data to match our binning by weighting 
according to the uncertainties in the LOSVDs. 

We re-construct Gauss-Hermite moments from the 
combined SAURON LOSVDs for plotting purposes only. 
Figure 13] shows these moments near the major and minor 
axes. The major axis V for the SAURON data is signif- 
icantly lower than that measured for the Gemini data. 
The reason for this is that SAURON data are binned 
to match the gridding of our model bins. Near the ma- 
jor axis, the bins range in polar angle from 6 = 0— 11°. 
These bins are described by a single LOSVD constructed 
by averaging individual LOSVDs which sample the re- 
gion at smaller spatial scales. Thus, the average LOSVD 
contains contributions from LOSVDs as much as 6* = 11° 
above the major axis. 

2.2.4. Globular Cluster Kinematics 



At large radii, we use individual glo bular cluster veloci- 
ties published in Bridg es et all (|2007| ) to derive LOSVDs. 
The data contain positions and radial velocities for 108 
globular clusters in NGC 4594. We discard the innermost 

14 GCs as there are too few GCs inside R < 130" to re- 
construct an LOSVD in the inner parts of the galaxy. As- 
suming axisymmetry, the positions of the GCs are folded 
about the minor and major axes. In order to preserve 
rotation, we flip the sign of the velocity for all GCs that 
are folded about the minor axis. The GCs are then di- 
vided into annular bins extending from ^ = 0° to 90° at 
radii of 131", 214", 350", 574", and 941" with roughly 
20 GCs per bin. 

Within each spatial bin, we calculate the LOSVD from 
the discrete GC velocities by using an a daptive kernel 
density es timate adapted f rom [Silvermanl (4986 ) and ex- 
plained in lGebhardt et al.l ()1996[ ). Each LOSVD contains 

15 velocity bins. The velocity bins are highly correlated 
for the GCs, and there are likely only a few degrees 
of freedom per LOSVD. The 1-cr uncertainties in the 
LOSVDs ar e estimated through bootstrap resamplings 
of th e data (jGebhardt et al.lll996l: iGebhardt fc Thomas! 
[20091) . 

We compute Gauss-Hermite moments from the GC 
LOSVDs — again for plotting purposes only — and show 
these in Figure [3] Uncertainties are calculated by fit- 
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Fig. 3. — Gauss-Hcrmite moments for NGC 4594 from various 
sources. Black diamonds with error bars are from Gemini long- 
slit observations along the major axis. Red diamonds are from 
SAURON data near the major axis. Light blue triangles are 
SAURON data near the minor axis. Green crosses are the 3 HST 
data points, and dark blue squares are from the globular clusters. 
Solid lines are the result of our best fit model. 



ting monients to each resampling of the LOSVD during 
the bootstrap. The GC kinematics resemble the minor 
axis stellar kinematics in many of the panels. For exam- 
ple, their dispersions appear to be an extrapolation of 
the minor axis velocity dispersions. There is slight rota- 
tion (small /13) and possible evidence of radial anisotropy 
(positive /14). 

3. DYNAMICAL MODELS 

The dynamical models rely o n the orbit superpos i- 
tion technique first developed bv iSchwarzschildl ()1979D . 
We assume axisymmetry and match the luminosity den- 
sity profile and kinematics of the galaxy to those recon- 
structed from an orbit library. The library is populated 



with orbits carefully chosen to sample E, L^, and the 
third, non-classical integral I3. 
The code used in this pa per is descri bed in 

(Thomas et all (pOOl 

Similar axisym- 

IRixetal.1 (11997') . 



Gebhardt et al 
20051) ^ 



,-- ''200(1 m. . 

and ISiopis et all (|2009H . 



metric codes are presented 



prese ntea 
] (fl998l) . |Ci 



van der Marel et all (|1998D . ICretton et all (119991) . and 
Valluriet al.l()2004[ ). Van den Bosch et al. ([2008!) present 



a fully triaxial Schwarzschild code. The basic outline of 
our code is as follows; (i.) convert the luminosity den- 
sity distribution ^{r) into the stellar density p{r) via an 
assumed stellar mass-to-light ratio M/Li. [ii.) Add to 
this density the contribution from a black hole of mass 
M, and a dark matter halo with density profile pDM{r). 
{Hi.) Calculate the potential $ associated with this den- 
sity distribution and integrate a large number of orbits 
(typically ~ 20,000) over many dynamical times, (iv.) 
Assign a weight Wi to each orbit and determine the Wi 
values by minimizing the x^ difference between the ob- 
served kinematics and luminosity density of the galaxy 
and those resulting from the PSF-convolved orbit library, 
subject also to the constraint of maximum entropy. 

We maximize the entropy-like quantity S = S — ax^ 
where S is the Boltzmann entropy and a controls the 
relative weight of S ov x^- For small values of a, repro- 
ducing the observed kinematics becomes unimportant, 
and the models act to only maximize entropy. As a in- 
creases, maximizing entropy becomes less important, and 
more weight is given to matching the observations. In 
practice, we start with a small value of a and gradually 
increase it until y^ a symp tote s. The interested reader 
may see ISiopis et all (|2009f ) or lShen fc Gebhardt! (|2010l ) 
for more details. 

Our model grid consists of 19 radial and 5 azimuthal 
bins covering a radial range of O'.'OS to 1800"spaced log- 
arithmically. Additionally, we use 15 velocity bins to de- 
scribe our LOSVDs. We incorporate the effects of seeing 
by convolving the light distribution for e ach orbit with a 
model PSF before comparing with data (Geb hardt et al.l 
l2000f l. We approximate the PSF as Gaussian with a 
FWHM of either 0'.'94, O'.'S or 0'.'09 depending on whether 
the data are from SAURON, Gemini, or HST observa- 
tions respectively. The convolution extends to a radius 
of 10 X FWHM. 

We run over 8,500 models with different values of the 
model parameters M/Lj, M,, and poM- We use Ax^ 
statistics to determine the best fit parameter values and 
their uncertainties. Models whose values of x^ are within 
Ax^ = 1 of the minimum for a given model parameter 
(marginalized over the others) define the l-cr or 68% con- 
fidence band of that parameter. 

3.1. Model Assumptions 

Our fiducial density profile is a combination of stellar 
mass, dark matter, and a central SMBH: 



p{r,e) 



M 

T 



y{r,e) + pDM{r,e) + M,5{r) (2) 



where M/Li is the stellar mass-to-light ratio, assumed 
constant with radius and 5{r) is the Dirac delta func- 
tion. The angle 9 is the angle above the major axis. 
While pdm can in principle be a function of 9, we do 
not consider flattened models. We assume a spherically 
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Fig. 4. — x^ ^'S a function of the four modeled parameters — M/Lj, Mm, Vc, and re- Every dot represents a single model. The solid line 
is a smoothed fit to the minimum, which represents the marginalized values. 



symmetric, logarithmic halo of the form: 



PDAiir) 



V' irt 



AtiG{7 



r2)2 



(3) 



This profile is cored for radii r <rc and produces a flat 
rotation curve with circular speed Vc for r ^ Tc- It has 
two free parameters, r^ and Vc, which are varied in the 
fitting process. Including M/Lj and A/,, this brings the 
total numb er of model parameters to fo ur. 

Recently, iGebhardt fc Thomas! (|2009[ ) have shown that 
the inclusion of a dark matter halo can significantly affect 
modeled BH masses. To test for this, we run a smaller 
suite of models without a dark halo. 

4. RESULTS 
Our best-fit values for the four model parameters are 



Mn 



M/L/ = 3.4 ± 0.05 f|^, M. 



376 ± 12 km s , and r 



1.6±0.4)x 10^ Mq, Vc 
4.7±0.6kpc. Figure H] shows 



the x^ minima around each of the model parameters. 
Each dot represents a single model and the solid curve is 
a smoothed fit to the minimum. The points of the solid 
curve at Ax2 = 1 above the minimum determine the 
1-CT confidence limits on the parameters. All four model 
parameters have well-behaved x^ curves with sharp, well- 
defined minima. This allows robust determination of the 
model parameters with small l-cr uncertainties. 

The uncertainties we present are derived strictly from 
Ax^ statistics. The l-cr error bars on quoted parame- 
ters correspond to models within Ax^ = 1 of the min- 
imum value. Systematic effects are likely to contribute 
in addition to this quoted uncertainty. While in gen- 
eral for other galaxies one of the biggest sources of sys- 
tematic uncertainty is inclination, for NGC 4594 inclina- 
tion uncertainties are unimportant. Other sources of un- 
certainty may include effects due to non-axisymmetries, 
but these are likely small or zero since only the 
most ma ssive ell i ptical s are thought to be significantly 
triaxial (|Binnevl Il978l: iKormendv fc lUingworthl 119821: 
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Fig. 5. — Correlation plots among the 4 parameters. Each dot represents a single model. Red dots are within the 95% confidence band 
and large black dots are within the 68% confidence band for an individual parameter 



iTremblav &: MerrittI I1996I ). For more on systematic 
uncert ainti es, the reader is referr ed to iGebhardt et alJ 
((200l and lGiiltekin"eFan (|2009allbl ). 

Figure [5] shows correlations among the four model pa- 
rameters. Plotted are the different projections of the 
four dimensional parameter space; every small dot cor- 
responds to a model run. Red dots are models that lie 
within Ax^ = 4 of the minimum, and large black dots are 
within Ax^ = 1- There appears to be a slight correlation 
between M/Lj and Af, — much less severe than in M87 
(JGcbhardt & Thomas 2009). Not surprisingly, the high 
resolution of our HST kinematics is able to break the de- 
generacy between M, and M/Lj. We discuss this further 
below. The dark halo parameters do not show any obvi- 
ous correlation, indicating the GC and stellar kinematics 
were able to break the degeneracy usually observed be- 
tween these two parameters. 

Our best fit model has (unreduced) x^ = 582.6. It 
is non-trivial to calculate the number of degrees of free- 
dom vdof- Roughly, i^DOF = Nlosvd X Nbin, how- 
ever there are complicated correlations between velocity 



bins (jGebhardt et alJ 120031 ). With this crude estimate 
for udof, our best fit model has reduced x^ — 0.6. 

We compare the modeled value of our stellar mass- 
to-light ratio with that obtained from evolutionar y pop- 
ulation synthesis models (Marastonj 11998 . 2005) . We 
adopt values of 10 Gyr a nd 0.1 for the stellar age and 
metallicity of NGC 4594 ()Sanchez-Blazauez et aITl2006D 
and use these to derive the predicted I-band M/Lj from 
the Maraston models. For a Salpeter IMF with stellar 
masses drawn from the range 0.1 — IOOMq, this analysis 
yields M/Li= 3.99. If instead the stars obey a Kroupa 
IMF drawn from the same range, then M/Li= 2.58. 
We multiply these by a factor of 1.096 corresponding 
to Aj = 0.099 t o correct for Galact ic extinction along 
the hnc of sight (jSchlegel et al.lll998[ ) to obtain M/Li= 
2.83 and M/ Li= 4.37. Our dynamically-derived stellar 
MjLi— 3.4 ± 0.05 falls nicely between these two values. 

In Figure [HI we plot the total mass-to-light ratio as a 
function of radius for our best fit model with 1-cr un- 
certainties (gray region). The red cross-hatched region 
represents the range in stellar M / Li from stellar popula- 
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Fig. 6. — Local dynamical mass-to-light ratio for the best fit 
model. The gray band indicates the 68% confidence band, as de- 
termined from the limits placed on the 4 model parameters. The 
red cross-hatched region indicates the extinction corrected stellar 
M/ Lj derived from population synthesis models. 



tion models described above. Total M /Li rises near the 
center of the galaxy due to the contribution of the su- 
permassive black hole. As we go out in radius, the stars 
become more important to the total mass over roughly 
the range 5" to 50". Here the total M/Lj approaches 
both our dynamically determined MjLi and the range 
derived from stellar population models. Past 50" MjLi 
once again rises due to the importance of the dark halo. 
Figure [7] plots the enclosed mass of each component as 
well as the total mass of the galaxy. At our innermost 
bin, the total mass is almost two orders of magnitude 
greater than the stellar mass, meaning we are probing 
the black hole's sphere of influe nce quite well. The green 
line p lotted is the mass profile iKormendv fc Westpfahll 
(|1989f) derive from their gas rotation curve. It agrees 
well with the total mass distribution derived here. 

4.1. Models without Dark Matter 

We run 189 models with no dark halo. In these mod- 
els, we exclude the globular cluster data and use only 
the stellar kinematics. Correspondingly, the number of 
degrees of freedom impacting the unreduced x^ ^^'^ pro- 
portionately fewer. We measure a black hole mass of 
M,= (6.6 ± 0.3) X 10* Mq and stehar mass-to-light ra- 
tio of M/Li = 3.7 ± 0.05. The minimum unreduced 
X^ — 628, proving models without a dark halo are a 
worse fit. 

We do not se e the dramatic change that 
IGebhardt fc Thomas! (|2009f) see in M87 where the 
inclusion of a DM halo causes their determination 
of M, to double. Ii istead , our results mirror those 
of iShen fc Gebha^ (|20Tol) in NGC 4649 where the 
inclusion of a DM halo does not significantly change the 
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Fig. 7. — Mass enclosed within spherical shells for our best fit 
model and 68% confidence region. The red line is the stellar mass 
profile while the black line and surrounding confidence region rep- 
resent the total mass (black hole -|- stars -|- DM). The dashed line 
is our bes t-fit dark matter halo. Green i ndicates the mass profile 
derived in IKormendv &: Westpfahll 119891 ) from gas rotation. 

modeled M,. The likely explanation for this behavior 
is the inclusion of high resolution HST kinematics in 
both NGC 4649 and NGC 4594. ISchulze fc Gebhardt! 
()2011|) find the same effect for a larger sample of 
galaxies. Whenever the data have high enough reso- 
lution to resolve the black hole's radius of infiuence 
i?inf ^ GMt/a"^, dark matter has no significant effect 
on the determination of M,. For NGC 4594 we measure 
i?inf ~ 57 pc ~ 1'.'2. We use HST/FOS kinematics 
whose central pointing has a PSF of 0'.'09 ~ 0.08i?inf. 
Additionally, the light profile of NGC 4594 is more 
centrally concentrated than that of M87. These factors 
combine to allow a more accurate determination of Af,, 
removing the freedom that the models have to trade 
mass between M, and M/Lj. This is evidenced by the 
lack of correlation among M, and M/Lj in Figure [5j 

4.2. Orbit Structure 

Having already determined the orbital weights that 
provide the best fit to the data, we reconstruct the in- 
ternal unprojected moments of the distribution function. 
We perform this analysis on our best fit model and the 
models that define the 68% confidence region (over all 
combinations of the four model parameters) , yielding in- 
ternal moments at each grid cell. 

We define the tangential velocity dispersion to be 

(Tt = \/ hi'^A ~^ '^e) where a^ is actually the second mo- 
ment, containing contributions from both streaming and 
random motion in the (j) direction. Figure [S] shows the 
radial run of the ratio ar/at- The second moment of 
the DF is tangentially biased where the disk is impor- 
tant (gray region) as expected but mostly is isotropic 
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Fig. 8. — Radial run of the ratio of the radial to tangential com- 
ponents of the velocity dispersion tensor. Shaded ares represent 
68% confidence regions, with gray indicating stars near the major 
axis, red meaning stars near the minor axis, and green representing 
GCs averaged over all angles. 

elsewhere. The red region plots Ur/ut for stars near the 
minor axis, showing almost perfect isotropy. The green 
region indicates that at large radii, globular cluster kine- 
matics show significant radial anisotropy. We discuss the 
implications of this in Section [5.21 below. 

5. DISCUSSION 

5.1. Black Hole-Bulge Correlations 

We discuss the position o f NGC 4594 on the M-a and 
M-L relat ions, as defined bv lGiiltekin et al.l (2009b, here- 
after [GO^ and compare our values of the correlation pa- 
rameters to previous measurements. We cal culate the 
effective velocity dispersion Ce similar to IG09I 



^l'{V\R) + a\R))I{R)dR 






.I{R)dR 



(4) 



where V{R) is the rotational velocity and I{R) is the 
surface brightness profile. This makes ae essentially the 
surface-brightness-weighted second moment. Instead of 
integ rating from the center of the galaxy {R = 0) as 
IG09I did, we integrate from R > Rinf to ensure that 
we do not bias o-g with the high dispersion near the 
black hole. Our outermost kinematic data point is at 
R — 45", thus we have a gap in kinematic coverage 
between 45" < R < Re = 114". The velocity dis- 
persion near the end of our long-slit data is dropping 
sharply, however V may still contribute to the integral 
for R > 45". To invest igate this, we use the ga s rota - 
tion curve presented in iKormendv fc Westpfahll ()1989f ) 
which extends well beyond Re- Truncating the integral 
Sit R = 45" gives ae — 292 km s~ while using the ex- 
tended rotation curve yields a^ = 297 km s~^. 

The problem with this definition of Ue is that it in- 
cludes a contribution from the rotation of the disk. It 



has been shown that black hole mass does no t correlate 
with disk properties ([Kormendv et al. 120111) so this is 
not ideal. However, to compare with IG09I we must be 
consistent in our calculation of a^- We therefore quote 
this value of a^, when we compare to the M-a relation 
determined by iG09l As we expect black hole mass to 
track bulge quantities, disk contribution to CTe is likely 
to add a source of intrinsic scatter to spiral galaxies in 
the M-a relation. In fact, spiral galaxies are observed to 
have larger scatter about M-a than ellipticals of similar 
ae- 

We also discuss some possible alternatives to ae where 
we attempt to remove the disk contribution. One option 
is to remove V"^ (R) from Equat ion Q altogether. Bulges 
are k nown to rotate, however, (jKormendv fc IllingworthI 
Il982t ) and this will likely underestimate ae- This crude 
calculation gives ae — 200 km s~ . 

Another option is to assume some degree of bulge ro- 
tation a_£riOTi^jf_we_assum^NGC4594 rotates isotropi- 
cally (jKormendv fc Illingworthlll982f) , then its flattening 
determines its pos i tion o n the V/a — e diagram (Binne^ 
Il978| ). iKormendvl (|1982f ) shows that the relation 



V 
a 



(5) 



approximates the isotropic rotator line to roughly 1% 
accuracy. We use our value of the bulge ellipticity e = 
0.25 in Equation ([5|) and assume this value of V/a applies 
globally to the entire bulge. We then use our measured 
dispersion profile a(R) to determine the bulge velocity 
H>uige(-R)- Using these quantities, we determine ae = 
23 km s~^. We compare this to the kinematics listed 
in IKormendv fc IllingworthI (|1982f) . These data include 
long-slit spectra taken at a position angle parallel to the 
major axis, but 30", 40", and 50" above it. From these 
data, it is apparent that the bulge a off the major axis is 
roughly constant at ~ 220 km s""'^. The rotation velocity 
rises from to 100 km s~ at large radii. We estimate 
the luminosity- weighted mean 1/ ^ 50 km s~^. Adding 
this in quadrature to the constant bulge a = 220 km s~ 
gives ae ~ 226 km s~^- This estimate does not contain 
any rotation from the disk, and is consistent with our 
determination of ae = 230 km s^ obtained by assuming 
a constant V/a- 

Our black hole mass Af.^ (6.6 ± 0.4) x 10^ M© agrees 
nicely with that of IG09I which uses the K88 value. In fact, 
when corrected for their various distance determina tions , 
most values of M, in the literature agree quite well ( K88l 
lEmsellem erallll994bl IK96I . iMagorrian et al.l 119981 ) de- 
spite the many modeling techniques and datasets used. 
This is likely due to the high degree of isotropy as evi- 
denced in Figure |8l deduced fr om the V/a — e dia gram 
(jKormendv fc iningwort'hlll982D . and noted in|K88|., 

Figure M plots the position of NGC 4594 on the IGQI 
Al-a and AI-L relations. We plot each determination of 
ae in the left-hand panel. Straightforward application 
of Equ ation (JSJ) leads to a value of ae that falls directly 
on the IG09I M-a line (blue asterisk). Next closest is the 
method of calculating ae by assuming a value of V/o' 
(green square). This point lies 0.44 dex above the IG09I 
line, however this is still within the estimated scatter. 
The calculation of ae that ignored all rotation is, not 
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Fig. 9. — Position of NGC 4594 on the IG09I M-cr and M-L relations. The plot of M-rr (left) shows the three ways we calculate CTe as 
well as the value from G09 (blac k tria ngle) . In order of increasing ere we plot Ce with no rotat ion (r ed diamond), ae assuming a value of 
V/a (green square) and ae as in |G09l (blue asterisk). For the M-L relation (right) we plot thelGOS value (black triangle) along with our 
measurement (green square). 

of tracer particles. This could also indicate the two pop- 
ulations have different formation scenarios. 

Unfortunately, there is no radius in the galaxy where 
we have simultaneous coverage of both stellar and GC 
kinematics. Thus, we are unable to test whether the 
stellar orbits become more radial in the ~ 50" between 
where the stellar kinematics run out and the GCs begin. 
However, since the light profiles of both populations are 
significantly different (Figure [1]) there is no reason to 
assume they should share similar orbit properties. 

Figure [5] shows the globular clusters in NGC 4594 are 
radially anisotropic {(Tr/at > 1) over roughly the ra- 
dial range 100-1000" (approximately 1-10 Re). Previous 
studies of the GC systems of gal axies have found their ve- 
locity ellipsoids to be isotropic ^ Cote et al. 2001, 2003). 
However, these studies used spherical Jeans modeling in- 
stead of the more general axisymmetric Schwarzschild 
c ode we use. 

IRhode fc Zepj (J2004I ) determine with high confidence 
that the color distribution of the GC system in NGC 4594 
is bimodal. This may indicate different subpopulations 
of GCs with different orbital properties that formed at 
different epochs in the galaxy's history. In our analysis, 
we make no distinction between red and blue subpopu- 
lations. We use the light profile and kinematics of all 
available GCs, regardless of color. However, since we 
use different sources for our kinematics and photometry 
data, there is the possibility that each source draws from 
a different GC subpopulation. 

5.3. Dark Halo 

The parameters Vc and Vc of our model dark halo 
imply a central dark matter density of pc = 0.35 ± 
O.IMq pc^'^ Using an improved Jeans modeling tech- 
nique, rRmpel fc Tenies (2006) model NGC 4594 and 
find a dark matter halo with central density pc = 
0.033 M0 pc~^, ten times lower than our value. They, 
however, measure a larger stellar M/Ljv = 7.1 ± 1.4 in 
the bulge. 

We plot the fraction of enclosed mass that is dark mat- 




' / Re 

Fig. 10. — Fraction M£,M/{Mi, + Mom) of enclosed mass that 
is dark matter as a function of radius. 



surprisingly, farthest from the IG09I line. Calculation of 
the relevant quantities for comparison with the M-L re- 
lation is straightforward, and we plot our v alue of M, 
and Ly (green square) along with that from IG09I in the 
right-hand panel. 

5.2. Globular Clusters 

As demonstrated in section [4.2) we find significant ra- 
dial anisotropy in the globular clusters. It is interest- 
ing that the stellar kinematics at smaller radii do not 
show this feature. This difference in orbital properties 
combined with the difference in their light profiles might 
suggest the GCs and stars are two distinct populations 
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ter as a function of half-light radius Re in Figure [TUl At 
1 Re there is already a roughly 50-50 mix of stars and 
DM. Inside of Re the dark matter still contributes a non- 
negligible fraction to the total mass content. 

In a study measuring dark matter properties in 1.7 x 
10^ local (z < 0.33) early- t ype g alaxies from the Sloan 
Digital Sky Survey, iGrillol ()2010t ) find a correlation be- 
tween the fraction of dark matter within Re and the log- 
arithmic value of Re- With our measured value of i?e, 
this correlation predicts a dark matter fraction at Re of 
0.68. Our value of 0.52 is smaller, but still within their 
6 8% confidence li mit. 

iThomas et al.l (|2009f) derive scaling relations for halo 
parameters based on observations of early-type galaxies 
in the Coma cluster. These relations are constructed 
for similar galaxies using the same halo parameteriza- 
tion and modeling code used in this paper. This makes 
comparison to our parameters straightforward. We com- 
pare to the observed relations between halo parameters 
Tc, Vc, and pc and total blue luminosity Lb- Our value 
of Vc falls directly on the Vc-Lb relation, however our 
measured re is smaller by roughly an order of magni- 
tude. Since pe oc Ve/re, the discrepancy in re causes 



our measurement of Pc to be high when compared to the 
I Thomas et al.l (|2009D Pe-Ls relation. Scatter in this rela- 
tion is large, however, and the environment of NGC 4594 
is different from that of the Coma ga laxies . 

iKormendv fc FreemanI ()2004l 12011^ also derive scaling 
laws for similar parameters in galaxies of later Hubble 
type (Sc-Im). We measure a m uch higher density and 
much s malle r core radius than the IKormendv fc FreemanI 
(|200l I20TI relations imply at the Lb of NGC 4594. 
We interpret this as the result of severe compres sion of 
the h alo by the gravity of the baryons (Blumentha l et al.l 
Il986( ). Such an effect is expected in early-type galaxies 
with massive bulges. 
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